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A semiclassical diagrammatic approach is constructed for calculating correlation functions of 
observables in open chaotic systems with time reversal symmetry. The results are expressed in terms 
of classical correlation functions involving Wigner representations of the observables. The formalism 
is used to explain a recent microwave experiment on the 4-disk problem, and to characterize the 
two-point function of the photodissociation cross section of complex molecules. 
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I. INTRODUCTION 



An experimental enquiry of the internal structure of 
a complex system involves, usually, some type of a scat- 
tering process. For instance, a photon is scattered from 
the system and then collected by a remote detector. In 
other situations, the collected objects are fragments of 
the initial system itself. A prototype example of the lat- 
ter process is the photodissociation of molecules: The 
molecule absorbs a photon, and disintegrate! by redis- 
tribution of energy in its vibrational modestJ. A com- 
mon feature of such systems is the coupling to contin- 
uum modes. Namely, above some energy threshold, the 
system is open. 

Open systems are characterized by resonances. These 
are eigenstates of open Hamiltonians which are normal- 
izable, and therefore correspond to complex eigenvalues 
associated with a decay in time. In most interesting sit- 
uations, this decay is sufficiently slow, allowing for the 
system to explore a large part of the phase space be- 
fore disintegration. In complex systems, such as nuclei 
or large molecules, the dynamics on these long lived res- 
onances is chaotic, therefore, the expectation value of a 
generic observable exhibits a statistical behavioitl. For 
instance, the absorption cross section for photodissoci- 
ation of large mokjcules is a pseudorandom function of 
the photon energytl'Q. This behavior suggests a statisti- 
cal analysis of observables in open chaotic systemda^llS. 

The main purpose of this paper is to construct a dia- 
grammatic scheme for calculating correlators of observ- 
ables in open chaotic systems. This diagrammatic ap- 
proach is similar to the cross diagram technic in disor- 
dered systemstj. However, although both techniques rely 
on the semiclassical approximation there are differences: 
The main advantage of the proposed scheme is in its ca- 
pability of describing individual systems rather than an 
ensemble of them. Ensemble averaging, as opposed to the 
energy averaging employed here, tends to erase "clean" 
features of individual systems. These features may have 
important manifestations, as will be demonstrated in this 
paper. On the other hand, disorder diagrammatics pro- 
vides a systematic way of calculating corrections due to 
quantum interference effects. At this stage we are not 



able to provide a similar prescription for general chaotic 
systems. Therefore, we confine our attention to open 
systems with decay time shorter than the time at which 
weak localization effects set in. The latter time scale, 
known as the EhernfesJ-time, diverges logarithmically in 
the semiclassical limitliil. 

The semiclassical analysis of a quantum system briiigs 
out the relation to its underlying classical dynamicsEj. 
For example, the classical dynamics of an electron in 
a disordered metal is diffusive, and ensemble averages 
of quantum observables of the electron are expressed 
in terms oL the spectral properties of the diffusion 
propagatorlij. In more general situations the classical 
evolution is described by the Perron- Frobenius pnerator 
whose spectrum, known as the Ruelle resonancesOEa, de- 
scribes the irreversible relaxation of probability densities 
in phase space. 

The classical spectrum of the system sets the impor- 
tant time scales of the problem. When the system is 
almost close there are two significant time scales: One 
is the decay time of the system, r^, which is inversely 
proportional to the typical width of the resonances. The 
closer is the system, the longer is the decay time. The 
second time scale, Tc, is the time which takes for a clas- 
sical density distribution to relax to the ergodic state on 
the energy shell, when the system is closed. In diffusive 
systems this time is known as the Thouless timdlj. 

A large separation between the classical time scales, 
Td S> Tc, implies a universal statistical behavior of the 
system on energy scales smaller than H/tc. Namely, the 
statistics is described by. a. random matrix theoryca suit- 
able for open systemsElrE3. However, as will be demon- 
strated in this paper, there are important manifestations 
of the nonuniversal behavior of the system, which are 
especially pronounced when Td is of the same order as 
Tc- The main use of the diagrammatic approach will be 
to calculate these individual imprints the system. This 
information is most important for constructing effective 
models from experimental data. 

The organization of this paper is as follows. In sec- 
tion 2 we derive the classical propagator of the system 
by energy averaging of Green functions. This propagator 
constitutes the basic building block of the diagrammatic 
approach developed in section 3. We shall restrict our 



considerations to systems with time reversal symmetry. 
In section 4 we present two applications of the formalism. 
One concerns tlie_microwave experiment on the iV-disk 
scattering systernc3. The second is the photodissociation 
absorption cross section of complex molecules. In sec- 
tion 5 we summarize our results, and mention directions 
of further studies. 



II. CLASSICAL PROPAGATION FROM 
QUANTUM GREEN FUNCTIONS 

The purpose of this section is to construct the building 
block of the diagrammatic scheme which will be devel- 
oped in this paper. It is the classical propagator, which 
is a generalization pOf the "diffuson" ladder diagrams of 
disordered systcmalj. Our construction will rely on the 
semiclassical apnEopdmation for the Green function of the 
quantum systemOo. 



A. The semiclassical Green function 

Let Ti be the Hamiltonian of an open system having d 
degrees of freedom, and 7Y(x) be its classical counterpart, 
where x — (r, p) is a point in the classical phase space. 
The advanced (+) and retarded (^) Green functions of 
the system are 



G±(6) = 



1 



e±iO-H 



where iO denote an infinitesimal positive imaginary part. 
In the semiclassical limit, these Green functions contain 
two contributions: 



G^ ~G 



w 



G 



± 

osc 



The first, known as the Weyl term, is a smooth function 
of the energy given by 



(r'|G±(6)|r> 



dp 



|-p-(r'- 



/i^ e ± iO - 7i:(x) ' 



(1) 



where h = 2irh is Planck's constant, and x = ((r -|- 
r')/2,p). The second contribution is an oscillatory func- 
tion of the jfinergy expressed as a sum over the classical 
trajectorieaHil from r to r' with energy e: 



(r|G_(e)|r') = E.^, 



*p--kSf,{r'.r-€) 



(2) 



Here Sfj,{r' , r; e) is the classical action of the /i-th trajec- 
tory, while A^ is an amplitude which can be expressed 
as a combination of second derivatives of S'^(r', r; e) with 
respect to r' and r. 



B. The classical propagator 

To construct the classical propagator of the system 
from an energy average of Green functions, consider the 
following function 

n(x',x;tj) = {Tr{G+ {e+nw)5{x-x)G- {e)S{x~x')}), 

(3) 

where the (5- function jOperator is defined by the symmet- 
ric Fourier transformEEl: 



5(x-x) = 



dqrfk^rk.(f_ 



l^2d 



^[k-(f-r)+q-(p-p)] 



(4) 



Although not written explicitly, it will be assumed that 
these (5-functions have a small finite width, say the inte- 
gration over k and q is limited to a large hypersphere. 
This width will be taken to zero at the end of the calcu- 
lation. Substituting definition (Q) in (g) we obtain 



n(x',x;c.) = |^e-i(PWV)^ 
G+(r' + |,r-|;, + MG-(r+|,r'-|;e: 



(5) 



Next, we substitute the semiclassical approximation for 
the Green functions. Each Green function has a smooth 
and an oscillatory contribution, thus, there are four terms 
in n(x', x; w). However, it is easy to see that two of them, 
{G^G'^g^) , vanish upon averaging. Therefore only the 
Weyl, {G+,G^), and the oscillatory, (G+^G^,^), contri- 
butions survive. 

The Weyl term is local in phase space, namely it is sig- 
nificant only when x ~ x'. This term, which we denote 
by H;oc(x', x; w), is of minor importance for our purposes, 
and we defer its calculation to Appendix A. The result, 
however, is 



H;oc(x',x;w) = 



2tt 



ih.h'^u 



-(5(x' - x),5(e - 7^(x)), (6) 



where uj~^ = uj -{- iO. 

Consider now the contribution from the oscillatory 
parts of the Green functions. The product, G+j^Gq^^, 
consists in a double sum over trajectories: /i from r— q/2 
to r'+qV2 and ly from r-|-q/2 to r'-q72, thus 



Gosc(.'^ + f^)Gosci'^) 



7 .^ti^l 



ASu,—iSiy 



(7) 



An illustration of two such orbits is depicted in Fig. 1. 
This term is, clearly, nonlocal in space. Its oscillatory 
dependence on the energy implies that upon energy av- 
eraging the main contribution comes from the diagonal 
part (/i = 1^) of the double sum ([7|) . The approximation 
of a double suroJpy a single sum is called the "diagonal 
approximation"Lj. It does not take into account quan- 
tum interference corrections such as weak localization. 
However, according to our assumptions on the relation 



between the Ehernfest time and the decay time of the 
system, these corrections can be neglected. 

E+h(D 




FIG. 1. An illustration of the classical trajectories con- 
tributing to the advanced (+) and the retarded (— ) Green 
functions in Eq. (0) before the energy averaging. 

The construction of the diagonal sum of orbits involves 
an expansion of the actions around that of the orbit going 
from r to r' namely: 



S,{r'- 



1 , 



flLOt,, 



Pp • q + Pu • q') 



^.(r'- 



^+|;e)~5.(r',r;e) 
1 , 

- o (p-- ■q + 



where p^ and p' denote the initial (at point r) and the fi- 
nal (at point r') momenta associated with the /x-th orbit. 
Similarly Pi, and p[, are the initial and the final momenta 
associated with the i/-th orbit, t^ is the time which takes 
for the particle to travel along the /i-th trajectory. Using 
the above expansions for the actions in (0), the diagonal 
approximation yields 

(G_(6 + nc.)G+,(6)) = Y. \A,\^e^^^--^+^'^-^+'^-'-). 

and integrating over q 

5{V-V,)5{V'-V',)- (8) 



Substituting this result into 
and q' we arrive at 



n(x',x;tj)= ^ 

orbits: r^r' 



A^\2^^ut^ 



This formula expresses n(x', x; lj) as a sum over the clas- 
sical trajectories from r to r'. The momenta (5-functions 
restrict the initial and final momenta of the trajectories 
to p and p', respectively. 

The diagonal sum ^ can be calculated using the flow- 
ing sum ruleE3 (proved in Appendix B): 



,d-lt;2 



- E 

orbits: r^r 



^Ml'5(x;„x^,tJ= (9) 

/•OO 

/ dtdp'dpg{x', X, t)S{e - H{x))6{x' - x(i)). 
Jo 



Here ^(x', x, t) is a general function of the space points x' 
and X, and the time t. The coordinates x^ = (r, p^) and 
x' = (r', p' ) are the initial and final phase-space points 
of the /i-th trajectory, while t^ is the corresponding pe- 
riod of the trajectory. Finally, x(t) is the phase space 
trajectory of the system, as function of the time t, start- 
ing from X. This trajectory is the solution of Hamilton's 
equations: r — 97i(x)/9p, p = ~dH{x.)/dr, with initial 
conditions x(t = 0) = x. Using (||) to replace the sum 
over orbits in (M) by an integral over the time we finally 
obtain 

f°°dt e*'^* 
H(x',x;..) = <5(6-7^(x))y^ __^(x'_x(t)). (10) 

The function n(x' , x; w) is the Fourier transform of the 
classical propagator of the system for t > 0. It is the gen- 
eralization of the "diffuson" of disordered systems to gen- 
eral chaotic systems. This analogy becomes clearer when 
projecting n(x', x; uj) down to real space, i.e. integrating 
over p and p'. The integration over p' is straightforward, 
and to integrate over p we assume that 7i(x) = p'^ /2m, 
and that r{t) is independent of p. This is good approx- 
imation for a diffusive motion on long time scales, since 
momentum relaxation is fast. Thus 



(ip'(ipH(x',x; 



2ttD 



dte"^'S{r' -r{t)), 



where D is the average density of states per unit volume. 
The function, 6{r' — r(i)), is the matrix element of the 
diffusion propagator between the position states |r) and 



S{r'-r{t)) 



-,DV' 



*|r), 



where D is the diffusion constant. Inserting a complete 
set of momentum states, which diagonalizes the propa- 
gator, and integrating over t we arrive at formula for the 
diffuson of disordered systems: 



/ 



(ip'(ipH(x',x;w) ~ -jy^}! 



„iq-(r'-r) 



Vn '^ -iuj + Da 
q 



2 ' 



(11) 



where V is the volume of the system. 

A spectral decomposition similar to (|ll| ) exists al 
the classical propagator of general chaotic systemsP"^ 
In the time domain one formally has: 



J(x'-x(i)) = (x'|e-^*|x) 

= (5(H(x) - H(x')) 



A^' 



(12) 
"^°*xL(x)x^(x'), 



where C = { ■ ,?i(x)} is the Poisson bracket opera- 
tor of the classical system, 7^ are its eigenvalues, and 
Xq(x) and Xq(x') are the corresponding left and right 



eigenf unctions. The above spectral decomposition re- 
sults from regularization of the classical dynamics. For 
instance, by keeping a finite width for the i5-functions in 
( [lO| ) and taking it to zero only after the diagonalization of 
the propagator. The resulting eigenvalues, called Ruelle 
resonances, constitute the Perron-Frobenius spectrum of 
the classical system. All classical eigenvalues of open 
systems have positive real part. They appear either as 
purely real or in complex conjugate pairs. We denote this 
set of eigenvalues by {7a}^0' ^^'^ order them according 
to the magnitude of their real parts, 7o < 7^ < 72 < • ■ ■, 
where 7a = 7a ±«7a- 

In concluding this section we remark that complex 
eigenvalues of the classical propagator, 7a = 7^ ± 17^, 
characterize ballistic chaotic systems. They emerge, for 
example, when typical classical trajectories spend long 
time near some short periodic orbit. This type of behav- 
ior does not appear in diffusive systems where classical 
relaxation is dominated by purely real eigenvalues. In 
section 4 we demonstrate the manifestations of complex 
Ruelle resonances in two examples. 



III. SEMICLASSICAL DIAGRAMMATICS FOR 
OPEN CHAOTIC SYSTEMS 



In this section we construct a diagrammatic scheme for 
calculating correlation functions such as the two-point 
function (Tr{G+(e -|- tiw)A}TT{G-{e)B}), where A and 
B are some observables, and the averaging, (••■), is over 
the "center of mass" energy e. The energy interval over 
which this averaging takes place is sufficiently wide to 
contain a large number of resonances, but narrow enough 
so that the classical dynamics within this interval is ap- 
proximately independent of the energy. 

The proposed diagrammatic scheme is a lift into phase 
space of the disordered diagrammatics which is embedded 
either in the momentum or real space. An important step 
in this direction is to express products of quantum ob- 
servables in terms of their Wigner representations which 
are functions of phase space variables. This issue will be 
considered in subsection A. Next we calculate the one- 
point (subsection B) , and the two-point (subsection C) 
functions. Finally (subsection D), the results will be gen- 
eralized to n-point functions by setting the diagram rules 
for their calculation. 



A. Wigner representations 



Let (rjAjr') be a matrix element of the observable A, 
where (r| and |r') are position states in Dirac notation. 
The Wigner representationcj. 



A(x) 



/ 



of the operator A is a function of the phase space vari- 
ables X — (r,p). It is a faithful representation of the 
quantum operator, since all its matrix elements can be 
reconstructed from the inverse relation 



(r|i|r') 



^,-^P(r'-r 



hd 



'A(x), 



(14) 



where x — ((r' +r)/2,p). Yet, expressed as a function 
on the classical phase space, the Wigner representation 
is convenient for semiclassical expansions. 

The external product of two operators A ® B, with 
matrix elements (ro|A|ri)(r2|-B|r3}, have more than one 
Wigner representation. These representations corre- 
spond to the various ways by which pairs of coordinates 
are used for the Fourier transforms. One possibility is 
to pair the coordinates of each one of the operators sep- 
arately, i.e. ro with ri and r2 with r3. Obviously, the 
result in this case will be the product of the Wigner rep- 
resentations oi A and B. Denoting it by [AB]d{x,x'), we 
have 

[AB]d(x,x') = A(x)B(x') ^ /"dqdq'e^(P-i+P'"i') 



x(r-||i|r+|)(r'-||B|r' + |), 



(15) 



where x = (r,p), and x' ~ (r',p'). The two other repre- 
sentations correspond to different pairing configurations: 
To with r3 and ri with r2, which leads to 



[AB]six,x')= /"dqdq'e^^Pi+P'"!') 

x(r-||i|r' + Hl)(r'-||S|r+|), 
and ro with r2, ri with r^, which gives 
[ylS]c(x,x')= /dqdq'e^^P-'i+P'-i') 



(16) 



(17) 



x(r-||i|r'-|)(r+||i?|r' + H.) 



dqe 



rpq/ 



';\A\v + ^), 



(13) 



Similar to the inverse relation (|lj), there are inverse 
relations for products of matrix elements. The simplest 
case is the inverse relation for [yli?]c((x, x') = A(x)_B(x'). 
It is a product of two inverse formulae as (M) . The in- 
verse relations corresponding to (Hq) and (|17D are: 

(ro|i|ri)(r2|S|r3) = 

= I ^[ABU^,^')e-^i^<^^'-'^+^'<^^'^-\ (19) 

where in M) x= ((r3-)-ro)/2, p) and x' = ((ri+r2)/2, p'), 
while in @ x= ((r2+ro)/2, p), and x' = ((ri+r3)/2, p'). 
In understanding the nature of the new Wigner rep- 
resentations, it will be instructive to consider examples. 
As a first example consider the operators 



^ = /,(f-ro), and B = /,(f-ri), (20) 

where /^(r) is a Gaussian function of width a: 



^'^^'■)=(2^°"p{"fe} 



(21) 



A straightforward calculation of the Wigner representa- 
tions yields 



[AB]d(x',x)=/,(r-ro)/,(r'-ri), 



(22) 



and 



[Ai?],^,(x',x) = h'^Siv - v')f,{v - ro)/,(r' - v^] 



x/r2 (p±p')exp 



(ri-ro) 



^^ +-(p±p')(ri-ro)|>, (23) 



where the + and — signs of p' correspond to [Ai?]c and 
[AB]s respectively. [Ai?]d(x',x) is independent of the 
momenta and large when r ~ ro and r' ^ ri. On the 
other hand, [A_B]s^c(x',x) are exponentially small for all 
values of x and x', when |ro - ri| 3> cr. If ro = ri and 
a|p|V?i2 » 1, then [AS],,,(x',x) « h'' fl{v - Vo)5{v' - 
r)5(p±p') 

As a second example consider the case 



i = S=|0)(0|. 



(24) 



Here both observables equal to the same projection op- 
erator. Assuming the wave function 0(r) to be real (as 
for the case of time reversal symmetry) one immediately 
sees that all types of Wigner transforms are precisely the 



[AA],,e(x',x) = [AA]d(x',x) - p^(x')/50(x), (25) 

where 



,(x)=y'dqe^P-'i(r-||</))(0|r + |) 



(26) 



is the Wigner function of the state (j){v). P0(x) is a 
real function which upon small smearing in phase space 
yields a positive definite functioruthat can be interpreted 
a "classical" density distributionCJ. 



B. One-point functions 

To begin with the calculation of correlation functions, 
it is instructive to consider, first, the simplest case of 
one-point functions: 

dx 



into the trace, and taking into account that energy av- 
eraging acts only on the Green function. This averaging 
leaves only the Weyl contribution ([|), therefore 



Ca=^((A(x))), 



(28) 



where double brackets denote the microcanonical average 
over the energy shell e = 7i(x), i.e. 

((A(x)))^ /'^A(x)5(6-H(x)). (29) 



Notice that A and ((yl(x)}) do not have the same dimen- 
sions, they differ by dimensions of energy. 



C. Two-point functions 

We turn, now, to the calculation of the two-point func- 
tion, 

CabH - (Tr{G+(e + ftw)i}Tr{G-(e)B})^, (30) 

where {'r^c denotes the connected part of the 
correlators. Within the semiclassical approximation it 
means that only the oscillatory parts of the Green func- 
tions contribute to Cab{^)- Inserting complete sets of 
position states one can write the correlator ( pO| ) as 

CABii^) = j dvodvidv2dY:iA{vQ,Vi)B{v2,V3) 

X (G+e(i"i, ro; e + ?iw)G^^^(r3, ra; e)), (31) 

where 

A(ro,ri) = (ro|i|ri), B(v2,v^) ^ {v2\B\v^). (32) 

Next, we evaluate the average of (G+ ^.G^^gj,). For this 
purpose notice that, to obtain a nonzero contribution, 
one has to pair orbits of the two Green function with 
similar actions. This condition imposes restrictions as 
for the possible configuration of the coordinates rp, ri, r2, 
and r3. The possibilities (in systems with time reversal 
symmetries) are: (i) ro ^ r^ and ri '^ r2, as illustrated 
in Fig. 2a; (ii) Tq ~ V2 and ri ~ rs, see Fig 2b; (iii) 
two possibilities corresponding to the case Tq ^ ri and 
r2 ~ r3, where the orbits are deformations of periodic 
orbits, see Figs. 2c and 2d. These pairing possibilities 
imply that the two-point function is a sum of terms: 



G. = (Tr{icjG-(e)})^y ^A(x)(cjG-(x)). (27) G..(.) = GS,(.) + G^,(.) + G-^(.), (33) 



Here A is an observable, and Q'G~(e) is the imaginary 
part of the retarded Green function. The second equal- 
ity in the above formula is obtained by substituting (llj) 



where the superscripts a, 6, c and d refer to the contri- 
butions associated with configurations of trajectories as 
shown in Fig. 2. 



(a) 





(c) 





(d) 



FIG. 2. The orbit configurations contributing to the aver- 
age of the Green functions product {G^sc(i'i,ro)G^3c(r3,r2)). 

We turn to calculate the various terms of Cab{'-^)- 
Consider, first, C'\q{lo) associated w\t\v the configura- 
tion of orbits shown Fig. 2a. We start by substituting 
the Wigner representation (|l^), implied by the pairing 
of initial and final points, into (pi|). Then we change 
variables to 



r = (ro+r3)/2, 
q = ra -ro, 



r' = (ri+r2)/2, 
q' = i"i -r2. 



and perform the average over the Green functions. 
This average gives precisely the classical propagator 
n(x',x;a;) (see Eq. H), therefore 



C1b{^) = /"dxdx'[^B],(x,x')n(x',x;^) 



(34) 



27r 



die-*(([^i3],(x,x(t)))). 



where in the second line we substitute formula ( px| ) for 
n(x',x;w), and integrate over x'. 

The calculation of the contribution, coming from or- 
bits of the type shown in Fig. 2b, follows along the 
same lines. The difference is that now we use the sec- 
ond Wigner representation, Eq. (|l^). Changing, again, 
variables to means and differences of the initial coordi- 
nates, and using the time reversal symmetry of the sys- 
tem, G'=^(r3,r2) = G=^(r2,r3), yields 

C%(w)= /"dxdx'[AB]c(x,x')n(x',x;w) (35) 

The computation of the last two contributions, associ- 
ated with the orbits of Figs. 2c and 2d, requires an addi- 
tional step. This step comes in order to impose the con- 
dition that the two orbits are deformations of the same 
periodic orbit. It is achieved by requiring that the orbit 
from ro to ri passes near (r2 -|- r3)/2, and similarly the 
orbit from ra to r2 passes in the vicinity of (ro-|-ri)/2. To 



implement this condition it will be convenient to intro- 
duce a local coordinate system r = (r, r j_ ) , where r is the 
time along the trajectory while the coordinate r^ is per- 
pendicular to it. Then the oscillatory part of the Green 
function satisfies the semiclassical product relation: 



G^ 



.(ri,ro;e) :^ h f dr±Gf^^{ri,r;e)G 



± 

osc 



(r,ro;e) (36) 



This relation can be proved by calculating the integral in 
the stationary phase approximation, see Appendix C. 

Substituting (^) in (^l]) and representing the matrix 
elements of the operators A and B in terms of their in- 
verse Wigner transforms ( [l4|) we obtain 



CTsi^) ^ J d^d^'[ABU^,x')K{^',^), 



(37) 



where 



/^(x',x) = — ^ /dqdq'drj^dr'Lff'e-«(P'i+P'i') x 

G+ c(r+|, r'+rl; e+nw)G+ir'+rl, r-|; e + aw) 



x = (r,p), and x' = (r',p'). In obtaining the above 
integral we also changed variables to 



r = (ro+ri)/2, 
q = ri -ro, 



r' = (r2+r3)/2, 
q' = rs - i"2. 



Now, we approximate the average of the four Green func- 
tions by the diagonal approximation. Diagonal sums in 
this case can be formed in two ways coming from the 
two possibilities of pairing the retarded and the advanced 
Green functions. These two possibilities correspond to 
the orbits configurations illustrated in Figs. 2c and 2d. 
The two contributions are identical due to the time re- 
versal symmetry of our system. As in the calculation of 
the classical propagator, we expand the actions to linear 
order around r, r', and e, and integrate over q, q', r^, and 



The result is 



i^(x',x) - 87r2/i2d^f'^ |A,j2|A,|2e^-(*. 



+*.) 



fiiy 



C/ P/J. ~r Pf N r-, 

X '5(p ^ )'5(Pm± - V>vi 



X ^(p' - 



P. 



)^(Pu± - P'lxJ 



where t^ and t 



are the periods of the orbits going from 
r to r' and from r' to r respectively, p^ and Pi, denote 
the momenta of the orbits at r, and similarly p' and p'j^ 
are the momenta at r'. We denote by subscript ± the 
perpendicular components of momenta that are conju- 
gate to r_L. Finally, we replace the sum over orbits by an 
integral using the sum rule (||), and obtain: 



Cfs^^) 



dtdt' 



e'-(*+*')(([AS]rf(x,x')))t,t'. (38) 



The average ((• • ■))t,t' of a general function, ^(x, x'), is 
defined as 

((5(x, x')))t,t' = j d^d^'gi^, x')<5(e-H(x))5(e-H(x')) 



X 5(xJ|-X||(i))J(x 



-xl|(i')), 



(39) 



where X|| denote a coordinate on the energy shell e — 
7i(x). The two (5- functions in the above integral imply 
that the average is over periodic orbits on the energy shell 
with period t + t'. The factor of two in ( |3^ ) is due to the 
time reversal symmetry. 

The final result for the two-point function, CAsi'-^), 
is the sum ( [3^ ) of (|4|), (p5|), and (pq). It expresses a 
quantum mechanical quantity, in terms of classical cor- 
relation functions associated with the Wigner represen- 
tations of the operators A and B. There are two sources 
of contributions to C^s(w): open trajectories, and peri- 
odic orbits. The magnitude of the corresponding terms 
depends on properties of the various Wigner representa- 
tions [AB]d(x,x'), [AB],{yi,x'), and [AB]c(x,x'). When 
A(x) and S(x') are very smooth functions, the main con- 
tribution comes from periodic orbits. If, on the other 
hand, A and B equal to the same projector, the main 
contribution comes from open orbits. 



D. Diagrams rules for n-point functions in systems 
with time reversal symmetry 

The results for the two-point function can be general- 
ized to n-point functions in systems with time reversal 
symmetry. It is instructive to formulate this generaliza- 
tion in terms of a set of rules for a diagrammatic calcula- 
tion. Only the connected part of the correlation functions 
will be considered here. It implies that the relevant part 
of the Green function is the oscillatory contribution (ph. 

We define the dimensionless propagator on the energy 
surface as 



n£(xj|,X|| 



= ^ / dt e 



Aut id 



/i'^<S(x||-Xi|(i)), (40) 



where both xii and xj, lie on the energy shell e = ?i(x). 
In order to shorten our notations, from now on, we omit 
the subscript ||. 

The rules for calculating the average of an n-point 
functions arc the following: 

• Write the correlation function as a space integral 
involving the matrix elements of the operators and 
the Green functions. 

• Find all possible pairs of initial and final points con- 
nected by Green functions, and construct the orbits 
configurations as in Fig. 2. To each coordinate pair 
assign a phase space point. 



• Express the Matrix elements of the operators as 
an inverse transform of the Wigner representation 
implied by the pairing configuration of the initial 
and final points. These Wigner representation are 
functions of the phases space points assigned in the 
previous step. 

• Phase space points connected by retarded and ad- 
vanced Green functions are associated with the 
classical propagator (|40|) connecting the two points. 

• In case where two classical propagators join at the 
same phase space point, then such a point is ac- 
companied by factor of l/27r (this factor emerges 
from expansions like (pq)). 

• Integrate over all phase space coordinates with the 
weight J(ix(5(e — ?i(x))(- ■ ■)/h'^- This integration is 
the same as microcanonical averaging (|29|). 

The diagrams of the two-point correlation function 
Ca_b(w) are shown in Fig. 3. The contributions C^^(w) 
and C\g{uj) are represented by Fig. 3a and Fig. 3b, re- 
spectively. These terms came from orbits configurations 
shown in Figs. 2a and 2b. The two other contributions 
(Figs. 2c and 2d), C^luj), are represented by the dia- 
gram of Fig. 3c. 

In Fig. 4 we show examples of diagrams contributing 
to the three-point correlation function. 



CABcii^i,iU2) = (Tr{G+(e + ftcc;i)i} 

X Ti-\^G-{e+nuj2)B\Ti(G-{e)c\\ 



(41) 



For instance, the top diagram of this figure equals to 
((ne(x", x', wi - W2)n,(x', X, u;i)[ABCUix, x', x")))/27r, 
where the microcanonical average is with respect to all 
the phase space coordinates. 




ng(x\x;co) 



(b) 



[AB] (x,x^) 



ng(x\x;co) 



[AB](x,x^) 




ng(x,x^;co) 



FIG. 3. The diagrams of the two-point function Cab{(^)- 
Wigglylines represent the classical propagator on the energy 
shell (kd), and dashed line correspond to the vertecies given 
by the appropriate Wigner representation of the observables 
A and B. 



+ e+ftcoi 



e x' e+fico2 




ne(x\x;(Bi) n£(x-\)f;C0,-C02) 

~^ lABCJ^UXx")'' 
nj(x",x\a)i-co2) 




nE(x\ X : ffli) ^"5^"^--/ n^Cx X ;Mi) 



FIG. 4. Examples of diagrams contributing to the 
three-point function ([fl|). The orbits configuration cor- 
responding to a diagram is drawn on its left side. 
[A_BC]a,^(x, x',x") denote various Wigner representations of 
the external product of observables, 4 (g) _B (g) C*. 

Finally we remark that, the procedure outlined above 
generates a large number of diagrams. However, knowing 
some general properties of the observables, such as sym- 
metries and smoothness of their Wigner representations, 
can help in reducing the number of diagrams consider- 
ably. 



IV. APPLICATIONS 



The purpose of this section is to present applications of 
the formalism developed above. In th n fixs t example we 
consider the iV-disk scattering systemEilE3. In the sec- 
ond application we study th&,indirect photodissociation 
process of complex moleculesEJ. Special attention will be 
given to understand the manifestations of the individual 
imprints of the systems. 



A. The Af-disk scattering system 



Consider a quantum particle moving in a two dimen- 
sional plain, and scattered from N disks located, say ran- 
domly near the origin, see Fig. 5. Resonances in this sys- 
tem can be pictured as situations where the particle is 
trapped for a long time in the vicinity of the disks. A 
natural quantity characterizing the system is the proba- 
bility of finding it in a final state 1-0/) when prepared in 
an initial state {ipi)- 




FIG. 5. An illustration of the classical TV-disk scattering 
system. In the experimental realization of the quantum ana- 
logue of this system, using microwave cavities, a transmitor 
located at ro excites a microwave field which is probed at ri . 



In the experimental realization of this systemtu, using 
microwaves cavities, a transmitor prepares the particle 
in a localized wave packet at rg, and a detector col- 
lects the particle at ri. The signal is proportional to 
|G+(ri, ro; e)p, and its simplest statistical characteristics 
is the correlation function: 



KH 



(|G'+(ri,ro;e)nG+(ri,ro;e + MP) 
(|G+(ri,ro;e)|2)2 



1. (42) 



To implement the diagrammatic approach, we start by 
writing |G''*'(ri, Tq; e)^ as a trace: 

|G+(ri,ro;e)|2=Tr{G+(e)<5(f-ro)G-(e)<5(f-ri)} . (43) 

The orbit configuration representing the above formula is 
shown in Fig. 6a. Assuming the distance |ri — ro|, to be 
larger than the particle wavelength there will be only one 
contribution associated with the oscillatory parts of the 
Green functions. The diagram representing this contri- 
bution is shown in Fig. 6b. A straightforward calculation 
of this diagram yields 

/■OO 

(|G+(ri,ro;e)P)= / dt /(n, ro;i), (44) 

Jo 

where 

/(ri, ro; t) ^ J J^S[e - H(xo)]<5[ri - r(xo, t)] (45) 

is proportional to the classical probability of finding par- 
ticle at ri, starting from Tq in all possible directions, and 
evolving according to the classical equations of motion. 
Here r(xo,i) denotes the trajectory of the particle, as 
function of time i, starting from the initial phase space 
point xo = (ro,po). 

The correlator K{lu) is proportional to the diagram 
shown in Fig. 6d, where again we assume the distance 
between ri and ro to be much longer than the particle 
wavelength. A straightforward calculation of this dia- 
gram leads to 



Kiu) 



F(ri,ro;uj) 



J^(ri,ro;0) 



(46) 



where 



F{lo) 



i^(ri,ro;cc>)= / di e'"*/(ri, ro; i) 
Jo 

is the Fourier transform of /(ri, Tq; i) for time i > 0. 
(a) (b) 



(47) 
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FIG. 6. The diagrams of the correlation function K(ijj). 
(a) The orbit configuration, and (b) the correspond- 
ing diagram for (|G'*'(ri,ro; e)| ) . (c) The orbit con- 
figuration, and (d) diagram for the connected part of 

(|G+(ri,ro;e)|'|G+(ri,ro;e-hftc^)P>. 

In understanding the typical behavior of K{lo), it is 
necessary to characterize the function /(ri,ro;<}. The 
spectral decomposition of classical propagator ( |l2|) im- 
plies thatE3 



/(ri,ro;i) ^ ^h^e 



-7q* 



(48) 



where 7q are the Ruelle resonances of the system, and 
ha are coefficients associated with the eigenfunctions of 
the classical propagator, Xa(x) and Xa('*^); ^s well as 
the precise positions of the transmitor and the detec- 
tor. The first eigenvalue 70 is the escape rate of the 
system. This eigenvalue dominates the behavior of the 
system at long times, i.e. small to. Taking only its contri- 
bution, the correlation function reduces to a Lorentzian, 
K{u)) K, 7o/(7o + '^^)n This result has been known long 
ago in nuclear physiesHj-Jt was reproduced in the context 
of chaotic scatteringMlSa, and-also measured experimen- 
tally in a microwave cavitiesEZI. 

We turn to consider the nonuniversal imprints of the 
system on K{lo). For this purpose we shall take into ac- 
count contributions of higher Ruelle resonances in (48). 
It will be convenient to rescale all frequencies with 70, 
and set 60 = 1 (without loss of generality) . Then, adding 
the next term of (Eq), F{ijj) takes the form 



1 + iuj 



± 



hi 



(49) 



where 7^ and 7" denote the real and the imaginary parts 
of 7i, respectively, and we assume 61 to be real. In Fig. 7 
we plot the correlation function K{uj) calculated from 
( p6| ) using (ES) . The solid line is the Lorentzian obtained 
when 5q = for a > 0. The dashed line is a representa- 
tive example of the nonuniversal behavior of the system 
(with 61 = 1.5, 7J = 2, 7" == 4, and b^ = for a > 1 ). 




10.0 



FIG. 7. The correlation function K{ul)). The solid line cor- 
respond to the universal limit where h^ = Q for all a > 0. 
The dotted line includes nonuniversal imprints of the system, 
with foi = 1.5, 7; = 2, and 7" = 4. 

A behavior very similar to that shown in Fig. 7 has 
been recently jObserved in a microwave experiment on the 
4-disk systemO. Our results explain the additional peak 
in the measured correlation function, -ft'(w), as the im- 
prints of the complex Ruelle resonances of the systemtfl. 



B. Statistics of photodissociation spectra 

As a second application of the semiclassical diagram- 
matic approach, we study the statistics of photodisso- 
ciation spectra of complex molecules, such as the rad- 
icals HO2 and N02- Disintegration of such molecules 
is a two step process. In the first step, a photon ex- 
cites the molecule to an energy above the dissociation 
threshold. Then fragmentation proceeds by redistribu- 
tion of energy in the vibrational degrees of freedom, or 
tunneling from binding to unbinding energy surfaces of 
the adiabatic electronic potentialEI. A barrier, which sep- 
arates quasi-stable states from continuum modes, hinders 
the immediate dissociation of the excited molecule. The 
large number of degrees of freedom and the complexity 
of the system imply that, on these long lived resonances. 



the dynamics of the system is chaotic. An illustration of 
such a system is depicted in Fig. 8. 




FIG. 8. An illustration for the electronic energy surfaces 
of methyl nitrite, CH3ONO, as an example for a complex 
molecule which undergoes an indirect photodissociation pro- 
cess. So and 5*1 are contour plots of the electronic surface po- 
tential associated with the ground and the excited electronic 
states. They are plotted as functions of the distances be- 
tween the Nitrogen and the Oxygen atoms, keeping all other 
degrees of freedom fixed at their equilibrium values. The ir- 
regular shape of Si leads to chaotic dynamics of the excited 
molecule. 

Consider a molecule, in the ground state \g), excited by 
a light pulse to an energy above the dissociation thresh- 
old, and let H denote the Hamiltonian of the system on 
the excited electronic surface (Si in the illustration of 
Fig. 8). It will be assumed that H represents an open 
system with several open channels. The photodissocia- 
tion cross section of the molecule, in the dipole approxi- 
mation, is given by 



a{e) =r;ImTr{iG-(e)} 



(50) 



where G (e) is the retarded Green function of the 
molecule, and A is a projection operator given by 



A = 10) (01, where 



T^\9)- 



(51) 



Here P = d • e is the projection of the electronic dipole 
moment operator of the molecule, d, on the polarization, 
e, of the absorbed light, and 77 = e/chsQ, c being the 
speed of light, and Sq the electric permitivity. The en- 
ergy e is measured from the ground state of the molecule. 
A natural statistical characteristic of the photodissoci- 
ation process is the dimensionless two-point correlation 
function: 



Ziu) 






It will be assumed that the excitation energy, e, is suf- 
ficiently high such that the mean spacing between the 



vibrational modes of the molecule is smaller than the en- 
ergy broadening due to the finite life time of the system 
in the excited states. This is the regime of overlapping 
resonances. 

From (pOl) , ([27|), and (pfl) it follows that the correlation 
function ( ^2|) equals to 



Z{uj) = if^—^2 — . 



(53) 



where A is the projection operator (p^), while Ca and 
Caa{i^) are the corresponding one-point and two-point 
functions. The special feature of the operator A ^ A 
is that all its Wigner representations are identical, see 
Eq. (Ea). It implies that (m) has the form 



Z{ijj) c± Zi{ijj) + Z2{uj), 



(54) 



where Zi{lu) is the contribution of open orbits (Figs. 2a 
and 2b), 



2 f°° 



((p^(x(0)p4x))) 
((P0(x)))2 ■ 



(55) 



while ^2(0;) is the contribution of periodic orbits (Figs. 2c 
and 2d), 



^2(0;) = 



^af? / dtdt e 

T^'^r? Jo 



f ^»c^(t+t') (W(x)p0(xO))t,t> ,gg- 



(W(x)))^ 



P0 (x) in the above formulae is the Wigner function of the 
initial state T>\g). 

Consider the limit of small cu. This limit reflects 
the long time behavior of the correlation functions 
((P0(x(i))p0(x))) and {{p^{x)p^{x')))t,t'- The leading 
contribution, in this case, comes from the smallest eigen- 
value in the spectral decomposition dl2), i.e. 70. Taking 
into account only this eigenvalue we obtain 

((p^(x(t))p^(x))) ^ ae-«'\ 

((p4x)p4x')))m' ^ 5e''^"(*+*'\ 

where 

a = ;i'^((p4x)x[,(x)))((p4x)xS(x))), (57) 

6 = ;i2'^((p4x)x[,(x)xS(x)))2. (58) 

The ratio b/a = A has an energy dimension. In the 
semiclassical limit, ?i — > 0, it approaches zero as h , thus, 
it corresponds to the mean spacing between resonances. 
Substituting the above approximations of the classical 
correlation functions into (p3) and (pq), and rescaling 
the energies as tiuj = fiA and hjo — TA, one obtains 



(52) zin)^z<^°\n) 



TT 



-t- 



1 



n^ 



T^ + n^ 27r(r2-h 172)2 



(59) 



where ^ = {a/{{P(/>)))'^ /b is a dimensionless constant of 
order unity. The first term in the brackets comes from 
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Zi{uj), i.e. from open trajectories of the type shown in 
Figs. 2a and 2b. This term constitutes the leading con- 
tribution to Z{uj) when F > 1. In the Hmit F ^ 1, 
the second term in the brackets becomes dominant. This 
term, associated with Z2{uj), comes from the periodic 
orbits of the system. Notice, however, that the regime 
F <C 1, corresponding to less than one open channel, is 
beyond our semiclassical approximation. 

Equation ( |59| ) represents the universal limit of the 
correlatcp|Z(ri). It was first derived rty Fyodorov and 
AlhassidEil using the nonlinear cr-modelE^. Our derivation 
confirms their conjecture that, in the limit of overlapping 
resonances (F > 1), Z{uj) can be derived by semiclassical 
methods. An alternative derivation of (|59|), using ran- 
dom matrix theory, is presented in Appendix D. 

Yet, the range of validity of formulae (|54| - p6|) goes far 
beyond the universal regime. They account also for sys- 
tem specific contributions which contain a valuable infor- 
mation about the nature of the system. This individual 
imprints come from the higher eigenvalues and eigenvec- 
tors of the classical propagator (^ . 

For example, the leading nonuniversal contribution 
coming from Zi{lu) is given by 



z[''m 



2^01 




^n 



(f7±f7i)2' 



10.0 



FIG. 9. The two-point function, Z{Q,), of the photodis- 
sociation cross section for systems with time reversal sym- 
metry. The solid line is the universal result (Eq. B9|) for the 
case 5 = r = 1. The dashed and the dotted lines show the 
two-point function which includes system specific contribu- 
tions in the approximation of Eq. (|60[). In both cases F = 1, 
Fi = 3 and ai = 6oi ~ 1, but, the dashed line corresponds to 
fii — 6, while the dotted line to f7i — 0. 



Here Fi = hj^/A and Qi = h^^/A, where j[ and 7" 
are the real and the imaginary parts of the second Ru- 
elle resonance. The constant Oi comes from an integral 
similar to ( p7| ) but with the eigenfunctions Xi''{^)- 
Similarly, the leading nonuniversal contribution of 

Z2in) is 



zi'\n) 



2Cfooi 



E 

± 



FFi -f7(0±f7i) 

(F2 + r!2)(F2 + (r!±r!i)2)^ 



where 601 is a constant depending on P0(x) and the eigen- 
functions Xo'^(''^) ^'^d Xi'^(x). 

The two-point function Z(ri), in the approximation 



z{n) ~ z(°)(fi) -I- z[^\n) + zi^\n), 



(60) 



is plotted in Fig. 9. The solid line represents the uni- 
versal form Z'^*'^(r2) given by Eq. (^9[). The dotted line 
shows the typical behavior of systems where decay of 
correlations is of diffusive nature. It appears when the 
subleading Ruelle resonances are purely real (fii =0). In 
this case the deviation from Z^'^' (f2) appears mainly as 
an increase of correlations near 17 = 0. The dashed line 
corresponds to the case of a complex Ruelle resonance 
(ill 7^ 0). It characterizes the behavior of ballistic sys- 
tems. Here the nonuniversal contribution is located in 
the tail of Z(rt), near fi = fii, where the universal term 
( p9| ) is already negligible. These plots demonstrate the 
significance of the individual imprints of the system on 

z{n). 



V. SUMMARY 

The main result of this paper is the construction of 
diagrammatic approach for calculating the n-point func- 
tions of open chaotic systems with time reversal sym- 
metry. In essence, this approach is simply an economic 
way for evaluating the "diagonal approximation" for the 
n-point functions. The result is expressed in terms of 
classical correlation functions associated with the observ- 
ables. These correlation functions reflect properties of 
the Perron-Frobenius spectral decomposition of the sys- 
tem, i.e. features of the irreversible dynamics of proba- 
bility densities on the energy shell. An important ingre- 
dient of this formalism is the generalization of Wigner 
transforms for external products of operators. 

The strength of this formalism is in its capability of 
treating particular systems rather than ensemble of them. 
Nevertheless, in the appropriate limit, our results reduce 
to those of random matrix theory. This was demon- 
strated by considering two examples: The A^-disk scatter- 
ing system, and the process of indirect photodissociation 
in complex molecules. The weakness of this diagram- 
matic approach is the lack of a clear systematic proce- 
dure for calculating quantum corrections such as weak 
localization correction. This limits our treatment to sys- 
tems with escape rate sufficiently large so that the system 
disintegrate before weak localization effects set in. 

The focus of this paper was on systems with time rever- 
sal symmetry, and no other discrete symmetries. Within 
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the universal limit, a generalization of our results to sys- 
tem without time reversal symmetry is straightforward. 
However, when considering nonuniversal manifestations, 
the situation is more complicated. The reason is that, 
usually in ballistic systems, there is no separation be- 
tween the time scales on which nonuniversal features ap- 
pear, and that where time reversal symmetry is broken. 
Further development of this diagrammatic approach 
should include a prescription for calculating weak lo- 
calization corrections, generalization to systems without 
time reversal symmetry, and to systems with other dis- 
crete symmetries, such as reflection or inversion. 
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Appendix A: The local term, n;oc(x',x; w)). 

In this appendix we calculate the Weyl contribution to 
the classical propagator, n(x',x;a;). Substituting ([|) in 
(^ and integrating over q and q' we obtain 

n;oc(x ,x;uj)= J ^^^ d{p )(5(p - p ) 



efi 



(po-Pl)-(r'-r) 



[e + hiu + iO- W(xo)][e ~ iO - n{xi)] ' 

where Xq = ((r -I- r')/2,po), and xi = ((r + r')/2, pi). 
Next, we change variables to k' = (po -I- Pi)/2, k — 
Po — Pi, and integrate over k'. The above integral, then, 
reduces to 



dk 



(5(p - p')e^''-('-'- 



h^'i [e + huj + iO- H(x+)][e -iO- H{x-)] ' 

where x^ = ((r + r')/2,p ± k/2). In the semiclassical 
limit, ?i ^ 0, the integral over k yields a 5- function, thus 



n;oc(x',x;a;) 



6{x' - :ii)/h'^ 



{e+hLu+io-n{x)){e-io-n{x)y 



(61) 



To extract the leading order semiclassical approximation 
of this expression, we manipulate the denominator as: 



den ■■ 



1 



huj + 


iO 


1 




huj + 


iO 


1 





1 



1 



e^iO-nix) 

1 



e + huj+iO-n{x) 
1 



e-iO-nix) e + iO-Hix) 

2mS{e-n{x)) 
nuj -\- i\j 

Substituting these results in (plh, we obtain m 



Appendix B: The classical sum rule (^ 

In this appendix we derive the sum rule (B). The 
derivation will be performed by calculating the right hand 
side of the equation and showing it equals to the left hand 
side. It is convenient to introduce a local coordinate sys- 
tem with time coordinate, r, along the trajectory, and rj^ 
perpendicular to the trajectory. The corresponding con- 
jugate momenta are the Hamiltonian function H ^ and 
p j^ . In these coordinates the right hand side of formula 
(P) takes the form: 



rhs — 



j dtj _^dp^dpY5(x, x', t)5{€ - H') 

X 5{H - H^)6{t - t[)5{v^ - v\t)5{vi_ - plJ. 

where r' and r are the absolute values of the velocity 
at the initial and final points respectively, while the sub- 
script t denotes the value of the corresponding coordinate 
after time t. Thus, in particular, r' ±t and p'j^^ are func- 
tions of the initial phase space point x' = (r', p') and the 
time t. Since energy is conserved, one has H[ = H'. The 
factor l/rr' in the above integral comes from the Jaco- 
bian of the transformation of variables. An integration 
over H, H' and px gives 



rhs = 



dt 



^5(r-r;)<5(r^-rlJ 
rr 



X 9{y^,y^,t)\H=H'=e, Pi=plj- 

Since the initial and final points of the particle as well 
as its energy are fixed, there is only a discrete set of tra- 
jectories which contribute to the above integral. These 
are the trajectories for which p'^ takes a value such that 
r^ = r'^j. We denote these trajectories by a subscript 
/i. Thus, p'^ and p^^ are the initial and the final mo- 
menta of the //-th trajectory, and i^ is the corresponding 
period. Then 



rhs 



M 



( dt 5{t-t^)f ^<5(ri-rl,)5(x^,x;,,i), 
Jo Jr., rr' 



where F^ is an infinitesimal region surrounding p'j_ in 
the momentum space, x^ = (r,p^), and xj^ = (r',pj^). 
The integral is straightforward, and the result is 



rhs = 2^g(x^,: 



,<p) 



TTyDet" 
rr 






where all quantities in the square brackets are calculated 
for the /i-th orbit. 

To see that this result is equal to the left hand side 
of Eq. (Q) we notice that the expression in the square 
brackets is related to the amplitude jA, of the semiclas- 
sical Green function (pi) according tccJ 



'^^'' " h^^^^ 



rr' \9p' ± 



(62) 
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Combining the last two equations yields the sum rule (|^). 

Appendix C: The product formula Eq. ( p6| ) 

In this appendix we prove the semiclassical formula 
( |36| ) relating a Green function to the product of two 
Green functions. We consider the case of G+g. The 
generalization for G~^^ is straightforward. We start by 
evaluating the integral on the right hand side of ( |M) us- 
ing stationary phase approximation. Substituting^) in 
(PI) yields 



rhs = nj: dr,rA,A,,e^ 



[S^(ri,r;e) + S^,(r,ro;e)] 



(63) 



where /i denotes the orbits from rg to r, while /i' the 
orbits from r to ri. The amplitudes A^ and A^i sat- 
isfy relation ( p^ . The stationary phase condition of the 
above integral is 

d 

0^ [S^,'{rl,r■,€) + Sf,{r,ro;e)] = 0. 

With the classical relations ^^^(r, rg; e)/9r = p^ and 
S'{ri,r;e)/dr = — p' , it implies that the momenta of 
the ii-ih and ^'-th trajectories at r are equal, p^ = p^/. 
Thus, r lies on a classical trajectory which goes from ro 
to ri. We denote this trajectory by i', see Fig. 10. No- 
tice that the stationary phase condition applies only for 
the perpendicular components of the momenta. However 
since the energy is fixed and equal for both trajectories, 
the longitudinal components of the momenta of both tra- 
jectories have to be equal as well. 
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FIG. 10. Classical paths joining ro ri and Tc -I- rx- 

Denoting by Tc a stationary value of r, we see that the 
stationary phase is the action of the z^-th orbit from rg 
to ri: 

5^'(ri,rc;e) -f 5^(rc, rg; e) = S'^(ri,rg; e). 

Thus, the integral m3t) takes the form 



r/is = ^/^e^^''('"^'''°''^\ 



where the pre-exponcnt factor, I,^, is given by 



I^ = fifAuA^th 



±lL (^ , fdp±^, p 



Dot 



9rj^ 



(64) 



It was calculated by expanding the sum of actions, 
S'p/(ri,r;e) -I- S'^(r,ro;e), to second order in r^ around 
re, and calculating the resulting Gaussian integral (|63|). 
To see the equivalence between this result and the left 
hand side of (EbI) it remains to show that I^, = Ay where 
Ay is the amplitude corresponding to the i^-th orbit from 
rg to ri . It is sufhcient to analyze the magnitudes of ly 
and Ay, since the Maslov phases are additive. Substitut- 
ing (62), for A^ and A^i , in (|6^) we obtain: 




Thus ly = Ay if the term in the curly brackets-equals to 
Det((9pj^i/9rj^g). Following Berry and Mounter we first 
write this term as a determinant of the product of the 
three matrices: 



{ } = Dot 



dY±_ 



dpi 



dpA 



9rj 



9rj 



(65) 



Now, we use the stationary phase condition P_l^ — P_l^' — 
0, which applies at r = re, and differentiate it with re- 
spect to r^i keeping r^g fixed. Then 



dpi 



9r_Li 



dr±^ 



9r^i 



Substituting this result in ( |65|) and using the chain rule 
of derivatives we obtain { } = Det(9p^g/9rj^i) which 
proves (P6|). 



Appendix D: An alternative derivation of ds 

In this appendix we derive formula ( p9| ) using random 
matrix theoryc^ (RMT). Our starting point is the Hamil- 
tonian of the form Tig -I- 17/2, where Tig is an A^ x iV ran- 
dom matrix of the GOE ensemble, and 7 is a constant 
which equals to the typical width of the resonances. With 
these definitions, a{e) — ((/)|(5.^/2(£ ~ 'Ho)\(j>), where 



i5^(a;) = — 



TT 7^ -|- x^ 



(66) 



and \(j)) — Vlg) (see Eq. ^l]). Denoting by faij) the j-th 
component of eigenfunction number a of Tig we have 






^(*=) = V2^2^S^/2{i ~ ia)^aij)^aik)(l){j)(t){k). (67) 
a jk 

To perform the ensemble averaging we use the property 
of RMT that eigenenergies Cq and the wave functions 
(faij) are statistically independent. Thus 
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He)) 



_n_ 

NA 



J2^ij)^ik)Sjk, 



into the expressions ( p7\ ) and ( p8| ) for a and b and evalu- 
ating ^ = a^/((p0))^/5, one obtains $ = 1. 



where 1/A — {J2a^i/'iie ~ ^a)) i^ ^^^ average density 
of states, and {ipa{j)^Pa(k)) = Sjk/N. Working in units 
where 'ri{(f)\(f)) — N, we obtain (cr) — 1/A. 

Consider now the connected part of the correlation 
function {a{e)a{e + uj)). Its calculation requires the av- 
erage of a product of two expressions like (^^ . Here we 
use the RMT relation: 

N^ {Va{j)Va{k)ipf3ii)(ppil)) = Sa(3[SjiSkl + 6jlSki] + SjkSil- 

The three terms on the right hand side of the above equa- 
tion correspond to the various paring possibilities of or- 
bits as shown in Fig. 2. In particular, the first and the 
second terms correspond to Figs. 2a and 2b, while the 
third term corresponds to Figs. 2c and 2d. Performing 
the sum over i, j, k, and /, the contribution of the first 
two terms is Iab{^) = 2{J2aS^/2ie-- (^a)Sj /2ie + uj -ea)), 
while the contribution of the third term is given by the 
correlator /crf(w) = {J2af3 ^'r/2{<^ - ea)Sj/2{e + uj- ep)) - 
1/A2. Thus, 



Iab{oj) 

and 



2 
A 



d€5^/2{e)6^l2{e + uj) 



Att 1^2 + 7^' 



Icd{^) = / d€de'5^/2{e)5^/2{e' - Lj)R{e - e'), 

where R{e — e') is the two-point function for the density 
of states of the GOE ensemble. Substituting the approx- 
imation i?(e — e') ~ K Jdiie**-^"^ ■'V^^j ^'^d performing 
the integrals over e and e' we immediately obtain 



2^2 ■ 



v.- JO 7r2 (72+^2) 

Collecting these results we arrive at 

Z(")(^) = A2(/,bM+/cd(c^)) 

7 



2A 7 A2 72-c^2 



TT Lj2 -|- 72 71-2 (724-1^2^ 

Formula (p^), apart from the prefactor ^, is obtained by 
rescaling the variables in the above equation as F = 7/A 
and fl = uj/A. 

We show, now, that in systems with a separation of 
time scales, 71 3> 70, the constant ^ in ( |59|) is approxi- 
mately unity. The separation of time scales implies that 
for times shorter than I/70 the dynamics of the system 
may be approximated by that of a close chaotic system 
whose Hamiltonian will be denoted by 7io(x). In this 
case the first eigenfunctions of the classical propagator 
are constant functions corresponding to the ergodic dis- 
tribution on the energy shell, namely x'(x) = x''(x) — 
[J dxS{e — 7io(x)]^^/^. Substituting these eigenfunctions 
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